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We develop and solve a continuum theory for the piezoelectric response of nanotubes under 
applied uniaxial and torsional stresses. We find that the piezoelectric response is controlled by 
the chiral angle, the aspect ratio, and two dimensionless parameters specifying the ratio of the 
strengths of the electrostatic and elastic energies. The model is solved in two limiting cases and the 
solutions are discussed. These systems are found to have several unexpected physical effects not 
seen in conventional bulk systems, including a strong stretch-twist coupling and the development of 
a significant bound charge density in addition to a surface charge density. The model is applied to 
estimate the piezoelectric response of a boron-nitride nanotube under uniform tensile stress. 
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I. INTRODUCTION 

Nanotubes and nanowires are a family of one dimen- 
sional structures whose properties can be tailored by 
their chemical composition, size and shape. The lat- 
eral dimensions of these structures can be reduced to the 
nanometer scale, a feature that is exploited in almost all 
applications envisioned for these objects. Carbon nan- 
otubes are the prototypical member of this family; they 
are the simplest chemically and are known to exhibit 
a diversity of electronic behaviors that are determined 
by the boundary conditions imposed on their electronic 
wavefunctions when they are mapped onto a cylindrical 
surface. While this sensitivity makes for rich theoretical 
study, it also makes CNTs difficult to use in applications 
since tubes must be identified and sorted according to 
their lattice structures. Boron-nitride nanotubes (BN- 
NTs) are the BI-V analog of CNTs and are an attractive 
alternative to CNTs because they have similar mechani- 
cal properties, are chemically inert over a wide range of 
conditions, and importantly, are all semiconducting with 
a band gap that is nearly independent of the chirality. 

A BN nanotube is a heteropolar structure with a lower 
lattice symmetry than the carbon nanotube. Impor- 
tantly, a twofold rotation about a perpendicular axis 
through the tube interchanges its A and B sublattices; 
this is a symmetry operation of the carbon nanotube but 
not of the boron nitride tube. In early theoretical work 
it was appreciated that this allows for a nonzero electric 
polarization in the ground state of the heteropolar tube 1 
and a linear coupling of the polarization to a uniform 
strain^, i.e. BN nanotubes are piezoelectric^. Interest- 
ingly, the piezoelectricity is a property inherited from 
the flat BN sheet, while the electric polarization of an 
unstrained tube is a new property that only occurs when 
the sheet is wrapped to form a nanotube. 

The description of piezoelectricity in a one dimensional 
nanotube differs fundamentally from the analogous for- 
mulation for a three dimensional crystal. In an ordinary 
three dimensional piezoelectric material described within 
the standard Landau-Devonshire theory, its macroscopic 



polarization (a vector) is linearly coupled to a uniform 
strain field (a tensor) through a third-rank piezoelec- 
tric tensor. Physically, this reflects the fact that bound 
charges arise from a divergence of the polarization, and 
these bound charges are mutually electrostatically cou- 
pled by the Coulomb interaction which diverges at small 
momentum proportional to 1/q 2 . Thus the free energy 
contains a term that is bilinear in the polarization that 
remains constant in the long wavelength q — > limit. As 
a consequence, a uniformly strained three dimensional 
piezoelectric induces a state of uniform bulk polariza- 
tion and the bound charges reside precisely at the sample 
boundaries. 

This state of affairs is upset for a one dimensional 
piezoelectric because the Coulomb potential that cou- 
ples the bound charges has only a softer — log q long 
wavelength divergence in one dimension. Thus the anal- 
ogous electrostatic coupling bilinear in the polarization 
vanishes in the long wavelength limit. The consequences 
of this are quite interesting and were partially explored 
in our previous worki In essence the strain and polar- 
ization must be treated as inhomogeneous fields, whose 
equilibrium forms can be obtained by minimization of 
an appropriate free energy that contains a nonlocal cou- 
pling between the one dimensional strain gradients. As 
a consequence, a uniform applied mechanical stress gen- 
erally induces a state of nonuniform strain in a nanotube 
or nanowire, and the bound charge redistributes from 
the tube end into its interior in a self consistent fashion. 
The distribution of the equilibrium strain and polariza- 
tion fields depends on a subtle interplay of their elastic 
and electrostatic interactions. Experimentally available 
nanotubes are expected to fall in a parameter regime 
where the elastic interactions dominate the problem so 
that most real materials, including BN nanotubes, are 
elastically "stiff." 

A second rather interesting consequence of this non- 
local theory occurs when different strain fields are sepa- 
rately coupled to the electric polarization. Elimination of 
the polarization then induces a bilinear coupling between 
different strain gradients. A striking consequence of the 
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induced coupling occurs on a nanotube where both tube 
extension (longitudinal strain) and tube torsion (trans- 
verse strain) couple to the polarization field. Here, a 
uniform tensile or compressive stress induces extension 
or compression and a concomitant torsion in its equilib- 
rium state. This conversion of a longitudinal to a tor- 
sion stress is impossible in a macroscopic homogeneous 
tube or cylinder, though it would generically occur in the 
elastic theory of a helical filament. We conclude that the 
appearance of an anomalous clastic response of this type 
is a macroscopic manifestation of the underlying (atomic 
scale) chirality of the structure. 

In this paper we develop a theory for this anomalous 
elastic response and apply it to study the equilibrium 
strained state of a family of BN nanotubes. We find 
that it has a rather rich structure, depending on the rela- 
tive strengths of the elastic and electrostatic interactions 
and on the chiral angle of the tube. Nonetheless, these 
anomalous couplings turn out to be fairly small, and will 
likely be difficult to observe directly in ordinary exper- 
imental environments. Physically, this reflects the fact 
that finite radius BN are mechanically stiff structures 
that tend to be dominated by their bare elastic interac- 
tions, weakly perturbed by the electrostatic effects of the 
type considered here. 

The rest of this paper is organized as follows. In SeclITl 
we develop a continuum theory appropriate for 1-D sys- 
tems and discuss limiting cases. In Sec. IIIII we present 
exact numerical and approximate analytic solutions for 
the various parameter regimes. In Sec. IIVI we apply our 
results to chiral BNNTs and calculate their piezoelectric 
response. Finally, in Sec. [V] we summarize our results 
and present conclusions. 



II. FREE ENERGY FORMALISM 

We consider a piezoelectric nanotube with radius R 
and length L, and let the tube axis be the z-axis. In 
what follows, we assume that L ^> R ^> ao, where ao is 
the B-N lattice spacing. In this limit of large aspect ratio 
and large tube radius a continuum theory is appropriate. 

The tube is placed under a constant external stress, 
which we partition into its longitudinal component f s 
and torsional component ft- The applied stresses in- 
duce strain fields r] s (z) (stretch) and n t (z) (twist). The 
induced polarization is related to these strains through 
P{z) = e s 7] s {z) + etrjt(z), where e s and e* are the linear 
piezoelectric constants that depend on the chiral angle 0. 
For a general chiral nanotube both e s and et are nonzero. 
For the two high symmetry classes of nanotubes only one 
piezoelectric coefficient is nonzero, e s for "zigzag" nan- 
otubes and e f for "armchair" nanotubes. The special 
cases where only one elastic degree of freedom couples 
to the polarization have no induced cross coupling terms 
and were considered previously.— 

The free energy of the system can be written as the 



sum of three terms, 

G = G c i as tic + Gdectro. _ 2_. I fi 7 li{ z )dz, (1) 

i J o 

where Gciastic is the elastic contribution, G e i ec tro. is the 
electrostatic contribution, and the last term represents 
the work done by the applied stress, and the sum runs 
over s and t. The elastic contribution is 

^ pL 

Gciastic = ^2 2 C ' lj J 1 li( Z ) 1 lo{ Z ) dz ' ( 2 ) 
i,3 

where the indices i and j run over s and t and we have 
introduced four 1-D elastic moduli Gy. The diagonal 
1-D elastic moduli are radius dependent quantities with 
units of energy/length, and are given by C ss = 2irRG\i 
and Ctt = 2ttRCqq. C\\ and Cqq are 2-D elastic constants 
in the conventional Voigt notation and are the relevant 
intensive variables, with units energy/area, which depend 
on the nanotube chemistry. The off-diagonal term, C s t, 
is small and decreases as 1/R 2 in the limit of large radius 
tubes. We therefore ignore this term by setting C s t = 
in all that follows. 

Using the fact that the bound charge density is re- 
lated to spatial derivatives of the strains through p{z) = 
— J2i e i9Vi( z )/9 z , the electrostatic contribution to the 
free energy can be written as 

G electro . =^ — JJ o —V(z-z )- Bp - dzdz , 

where again the sums on i and j are over s and t. Here, 
Viz — z') is a Coulomb kernel describing the electrostatic 
interaction between rings of charge centered at positions 
z and z' along the tube axis. e s and et are 1-D piezo- 
electric constants with units of charge. They may be 
obtained by rotating the piezoelectric constants for an 
infinite 2-D sheet onto the symmetry axes of a tube. 3 
The three-fold symmetry of an infinite hexagonal lattice 
reduces the number of independent 2-D piezoelectric con- 
stants to one, which we take as ei = e xxx with units of 
charge /length. If we define the chiral angle as the angle 
between the nanotube axis and a 2-D primitive trans- 
lation vector (6 = 30° is the zigzag configuration and 
9 = 0° is the armchair configuration), then the 1-D piezo- 
electric constants are given by 

e s — Ge2sin(30) 

e t = -Ce 2 cos(30) , (4) 

where G is the nanotube circumference^ These expres- 
sions are derived neglecting finite radius corrections due 
to tube curvature. However, these expressions agree with 
numerical results to within 15% for the smallest nan- 
otubes considered, with better agreement for larger ra- 
dius tubes where curvature effects are reduced. 
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The polarization induced along the tube is composed of 
dipoles with both a longitudinal and an azimuthal com- 
ponent. The longitudinal component accounts for the 
electrostatic energy given in Eq. ©. The total dipole 
moment of the azimuthal polarization integrates to zero, 
but higher order moments do not vanish and may con- 
tribute to the energy. The energy from these interacting 
rings of dipoles is of the form 

G r ing S ~ - J J dxdx'rji(x)K(x - x')rjj(x'), 

where the integral is written in terms of the dimensionless 
variable x — z/R and ef is the appropriate piezoelectric 
constant giving the local azimuthal polarization in terms 
of the induced strains. The important result is that the 
kernel goes asymptotically as K(x) ~ l/x 2m+1 , where 
to is approximately the number of 2-D unit cells in one 
tube circumference. This kernel decays so fast that we 
may approximate it with a delta function. This term 
then has the same functional form as Geiastic- However, 
as we will see below, the terms we will keep in G e i as ti c 



scale with the radius as R 2 , whereas all terms in G r i ngs 
scale as 1 /R 3 . Since we are interested in the limit of large 
radius tubes, this energy represents a minor correction to 
the elastic energy of the system and will be ignored. 

The strain field may be discontinuous at the bound- 
aries of the nanotube. To account for this, we let 
T]i(z) = gi(z)H(z), where gi(z) is a continuous function 
of z and H(z) = 1 for < z < L and otherwise. Writ- 
ing rji(z) in this form is convenient because it allows us to 
keep track of the delta function contributions to r][(z) at 
the boundary of the nanotube. In particular, it removes 
any ambiguity that might arise in an expression such as 
77^(0). In what follows we will call both rji(z) and gi(z) 
the strain field, since they agree over the length of the 
tube. 

We extract the dimensional dependence of the energies 
in Eq. ([1]) by expressing all lengths in units of the tube ra- 
dius and all radius dependent quantities by their relevant 
intensive quantities. We define a scaled length x = z/R 
and the dimensionless aspect ratio as Ai = L/R. Then, 
after integrating the electrostatic free energy by parts, 
the free energy is given by 



G = Tri? 2 f dx[G ll5s 2 (x) + G 66 .g t 2 (x)] -R [ dx[f s9s {x) + f t g t {x)] 
Jo Jo 

f Xl d 2 V(x-x') 

+2n 2 Re 2 2 / dxdx' v ' [sm 2 {36)g s (x)g s (x') - 2cos(30) sm(36)g s (x)g t (x') + cos 2 (39)g t (x)g t (x')] 

-An 2 Re 2 2 V'{Q+) / dx[sin 2 (3% 2 (x) - 2 sin(30) cos(38)g s (x)g t (x) + cos 2 (30) 5 2 (x)] , (6) 
Jo 
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where V'(0+) indicates the one-sided right derivative of 
the Coulomb kernel. It may be that V'(0+) = 00, but 
there is an equal and opposite infinity in the double in- 
tegral so that the energy always remains finite. 

The equilibrium strain fields are obtained by minimiz- 
ing G with respect to the strains and satisfy 

f Al d 2 V(x - x') 
g(x) - 2V'(0+)Ag(x) + J dx' ^ 1 Ag(x') = (5. 

(7) 

For notational simplicity we have defined a strain "vec- 
tor" as g{x) — (g s (x), g t (x)), and an elastic limit "vector" 
as (3 — (P s ,(3t), where (3 S — f s /(2nRCn) is the value of 
the uniaxial strain in the absence of piezoelectric effects 
and (3t = /t/(27ri?C66) is the value of the torsional strain 
in the absence of piezoelectric effects. The matrix A is 
given by 

/ A 2 sin 2 (3(9) -A 2 sin(30) cos(36») \ 
W \ -A3 sin(30) cos(30) A 3 cos 2 (30) J ' 

where we have introduced two dimensionless constants, 
A2 = e\/{RC\i) and A 3 = e^/{RC^), which measure 



the relative strengths of the electrostatic and elastic re- 
sponse. Note that C\x and Cqq are generally of the same 
order of magnitude, so A2 and A3 will also be of the same 
order of magnitude. 

We will see below that the solutions to Eq. ([7]) asymp- 
totically approach the elastic limit with a characteristic 
decay length controlled by A2 and A3 . There are two lim- 
iting cases of interest. When the tube is much longer than 
the decay length (Ai 3> 27rA2,27rA3), the strains are near 
their elastic limits over most of the tube length, with de- 
viations from the elastic limit near the tube ends. Such 
systems are dominated by their elastic response to the 
applied stress, and we call this the elastically dominated 
limit. Conversely, when the tube is much shorter than 
the decay length (Ai <§; 27rA2,27rA 3 ), the strains never de- 
cay to their elastic limits. Such systems are dominated 
by their electrostatic response to the applied stress, and 
we call this the electrostatically dominated limit. Solu- 
tions for these two limiting cases will be explored in the 
next section. 

The matrix A is not invertible, therefore the strains 
are not linearly independent. We find that the strains 
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are related by 

/x * A3 008(36*) . „ ... ,„. 

Plugging this relationship into equation (JT)) allows us to 
write the equilibrium equation in terms of g s (x) alone, 

g,(x) = A-2V'(0+)[79,(i)-t] 

where we have defined two new constants, 

7 = 27r(A 2 sin 2 (3<9) + A 3 cos 2 (36»)) (11) 

and 



r = 27rcos(36»)[A 3 cos(39)/3 s + A 2 sin(36»)/3 t ] 



(12) 



In order to solve Eq. (JTUJ) an appropriate Coulomb ker- 
nel must be specified. In the simplest continuum model 
where the kernel describes the interaction of charged 
rings at positions x and x' on a tube, the Coulomb kernel 
is given by 



V{x-x') 



1 

2^ 



d6 



c') 2 + sin 2 (6») + (1 -cos(6>)) 2 

This kernel diverges logarithmically cLS x — ^ X ' and leads 
to an infinite energy for surface charges. Using this kernel 
the boundary condition on the strain fields is e s ?y s (0) + 
e* 7 ?* (0) = 0. A divergent energy is an artifact of our 
continuum model and signals the breakdown of the model 
at short distances. It is useful to replace the divergent 
kernel with a softened kernel, 



V so i t (x 



1 



a 



(14) 



which retains the correct long range behavior but re- 
mains finite clS x — ^ X . This kernel allows the presence 
of a nonzero surface charge, with the "cost" of a sur- 
face charge controlled by the cutoff a. The tube radius 
forms a natural length scale below which the continuum 
model breaks down, suggesting a « 1 is an appropriate 
cutoff. In Rcf. |4j) it was shown that using the softened 
kernel with a = 1 gave nearly identical results for the 
strain, charge density, and potential across the tube as 
those obtained using the divergent kernel. The only no- 
table difference is that the softened kernel gives a non- 
zero surface charge which is replaced by a logarithmically 
divergent bound charge density at the boundaries with 
no surface charge for the divergent kernel. These are 
indistinguishable for all practical purposes. 

Once the charge density is computed, the potential 
along the tube is determined by 



U(z) = V(z- z')p{z')dz' 
Jo 



(15) 



After expressing all lengths in terms of the tube radius 
and integrating by parts, the potential becomes 



tt, \ e * f Xl , ,sdV{x-x') , , 



dx' 



(16) 



Differentiating this expression and writing all constants 
in terms of intensive variables allows us to relate the po- 
tential to the strain through 



dU(x) 
dx 



e 2 



A 2 sin(30) 



(17) 



Integrating this expression give the potential as an inte- 
gral over the strain field, 



U(x) 



(-2 



A 2 sin(30) 



Ai/2 



dx'(g s (x')-(3 s ) 



(18) 



where we have set the potential to zero at middle of the 
tube. 

Once gt(z) is obtained through Eq. ©, the change 
in the azimuthal angle, 4>, at a point z is given by 
Rdcf) = g t {z)dz. The total twist induced along the tube 
is therefore given by 



A</> = / dxg t (x) , 



(19) 



where we have transformed to the scaled length. Fol- 
lowing the usual convention, the total twist is posi- 
tive if the tube rotates in the counterclockwise direc- 
tion, and negative if it rotates in the clockwise direc- 
tion. Notice from Eq. (J9j) that for a general chiral tube 
there is a coupling between longitudinal and torsional 
strain such that a uniaxial stress alone ((3t = 0) is suf- 
ficient to induce twisting in the nanotube. This cou- 
pling vanishes in the high symmetry zigzag and arm- 
chair nanotubes, as may be seen by writing Eq. © as 
A 2 sm(39)(g t (x) - p t ) = A 3 cos(36>)(/3 s - g s (x)). For arm- 
chair nanotubes sin(36') = 0, which forces g s (x) = (3 S but 
sets no restrictions on gt(x). Likewise, for zigzag nan- 
otubes 008(36*) = 0, which forces gt(x) = (3t but sets no 
restrictions on g s (x) . For a general chiral tube the origin 
of this "stretch-twist" coupling lies in the electrostatic 
term in the free energy. In regions where one flavor of 
strain has large spatial gradients the bound charge den- 
sity due to that flavor of strain is also large. The overall 
bound charge density may be lowered in that region by 
gradients in the other flavor of strain that produce bound 
charges of opposite sign. The equilibrium strain fields 
are those that create a bound charge density profile that 
gives the lowest electrostatic free energy consistent with 
Eqs. ([51 and (TiT)|) , 

If we consider a system with (3 t = 0, then Eqs. ([9]) and 
(P| give 



= c 



cos(36>) 
sin(36>) ' 



(20) 
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where C = (A3/A2) J (/3 S — g s {x))dx. We will see below 
that for a system under tensile stress (i s > g s {x) for all 
x. In this case C > 0, and the sign of the total induced 
twist in the nanotube is entirely determined by the sign 
of the chiral angle (mod 27r/3). We find that A(j) > for 
right-handed nanotubes and A0 < for left-handed nan- 
otubes. Therefore, the sign of the induced twist may be 
used to determine the handedness of a chiral piezoelec- 
tric nanotube. There is no known order parameter for the 
handedness of a chiral nanotube that can be built up out 
of the two primitive vectors of the 2-d hexagonal lattice 
and the translational vector of the nanotube. As such, 
an unambiguous assignment of nanotube handedness is 
difficult with structural information alone. A measure- 
ment of the sign of the stretch-twist coupling allows one 
to easily make this assignment. 

A non-piezoelectric chiral nanotube will also couple 
stretch and twist to some degree due to the rearrange- 
ment of the atomic basis upon stretching. This effect is 
absent in our formalism because we have set C s t = 0, 
when in reality this quantity is small but non-zero for 
a chiral nanotube. It is worth noting that the electro- 
statically induced stretch-twist coupling described above 
is in general several orders of magnitude larger than the 
coupling in non-piezoelectric nanotubes. The torsional 
strain induced in a piezoelectric system can easily be of 
the same order of magnitude as the uniaxial strain. 



III. EQUILIBRIUM SOLUTIONS 

The equilibrium solutions for the strain fields gi(z) 
(Eqs. [S] and [TO]), the bound charge density p(z) = 
~J2i e i^ r li( z )/9z, and the electrostatic potential U(z) 
(Eq. I18p were determined using finite element analy- 
sis to solve Eq. ([TO)) . Convergence of the solution was 
monitored by increasing the number of finite elements in 
increments of 100 until the difference in successive solu- 
tions was less than 1% at the end of the tube where the 
convergence is the slowest. 



A. Elastically Dominated Limit 

The elastically dominated limit is defined by the con- 
dition that 27rA2,27rA3 <C \\. We investigate the case of 
a tube under constant uniaxial stress in the absence of 
torsional stress. The equilibrium strain fields for several 
chiral angles are presented in Fig. [TJ The strain fields 
rapidly relax to their elastic limits, and over much of 
the tube rj s {x) « (3 S and rj t (x) w 0. Near the bound- 
aries of the tube the electrostatic interactions effectively 
increase Cn, suppressing the longitudinal strain. Impor- 
tantly, notice that for chiral tubes r\ t (x) ^ for all x even 
though f3t = 0. Because the longitudinal and torsional 
strains couple through spatial derivatives the chiral tubes 
develop some nonzero twist near the tube ends. The total 




FIG. 1: Equilibrium solutions for a piezoelectric tube in the 
elastically dominated limit using the soft kernel. The upper 
panel shows the induced uniaxial strain for tubes of various 
chirality, while the lower panel shows the corresponding in- 
duced torsional strains. The graph is broken along the x-etxis 
in order to show details near the tube boundary. The graphs 
were generated using (3 S — l,/3 t = 0, Ai = 1000, A2 = IO/71-, 
A3 = 2A2,and e2 = y/E/n (charge/length). 



twist induced in the tube reaches a maximum at 9 = 15° 
and goes to zero as 9 — > 0° or 9 — ► 30°. 

The bound charge density and the electrostatic poten- 
tial are plotted in Fig. [2l The charge density is localized 
at the ends of the tube and decays rapidly on the inte- 
rior. Not shown are delta functions in the charge density 
right at the tube boundary. These are surface charges 
which appear because the polarization is nonzero at the 
tube boundary with the soft kernel. The potential varies 
rapidly near the tube ends due to the presence of the 
induced bound charges. On the interior of the tube the 
potential varies slowly and resembles the electrostatic po- 
tential produced by two opposite point charges at the 
ends of the tube. Almost the entire potential drop across 
the tube occurs near the tube boundaries. 

Both the bound charge density and the potential are 
maximized by letting 9 = 30° and decrease monotoni- 
cally as 9 — > 0. We find that the total potential differ- 
ence across the tube is well fit by AU ~ sin(30). There- 
fore, under longitudinal stress, the piezoelectric response 
is insensitive to small perturbations around the zigzag 
configuration. 

A striking feature of Figs. [1] and [2] is that the effects 
of the chirality only manifest themselves at the ends of 
the tube. The bulk of the tube responds to an applied 
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FIG. 2: Equilibrium solutions for a piezoelectric tube in the 
elastically dominated limit using the soft kernel. The up- 
per panel shows the induced bound charge density for various 
chiral angles and the lower panel shows corresponding elec- 
trostatic potential along the tube. The graph is broken along 
the i-axis in order to show details near the tube boundary. 
The graphs were generated using j3 s = l,/3 t = 0, Ai = 1000, 
A2 = 10/-7T, A3 = 2A2,and ej. — ^/E/tt (charge/length). 



stress as an ordinary elastic medium. This is a somewhat 
surprising result because it implies that boundaries are 
necessary for chiral effects to appear and that the bulk 
of the tube has no "knowledge" of its own chirality. Al- 
though this is evident from our free energy formalism an 
explanation of this effect based on the symmetry proper- 
ties of the tube is desirable and is the subject of current 
investigations. 

An approximate analytic solution to Eq. (jTUJ) may be 
obtained with a perturbation expansion and gives 



g s (x) =/3 a - (j/3 a - r) 



1 



1 



(x + a) 2 (Ai — x + a) 2 



This expression always gives the correct asymptotic de- 
cay to the elastic limit far from the tube ends. If 
the system is very far into the elastic limit, so that 
(2^X2/ a 2 ) <C 1, then this expression holds over the en- 
tire length of the tube. Using Eq. (fT8| to compute the 
potential difference, AU = U{X\) — {7(0), across a tube 
under uniaxial stress only gives 



AU = C L 



47re 2 sin(30)/3 s 
Ai + a 



(22) 



The constant Cjj can be determined numerically. If the 
system is extremely elastic {{211X2/ a 2 ) <C 1), then Cjj = 



47re2 sm{36)[3 s /a. 

We may use Eq. ^ to obtain gt(x) from g s {x), 
and then compute the total twist along the tube using 
Eq. (|19p . For a tube under uniaxial stress only, this gives 



A(f> = C, 



<t> 



4ttA 3 sin(30) cos(30)/3 s 



Ai 



a 



(23) 



The constant can be determined numerically. If the 
system is extremely elastic ((27rA 2 /a 2 ) -C 1), then — 
4ttA 3 sin(30) cos(30)/3 s /a. 

Equations (|22f and (|23|) show that, for systems in the 
elastic limit, the induced potential and the total twist 
essentially do not scale with tube length. This is a con- 
sequence of the fact that the bound charge density is 
localized near the tube ends with a decay length set by 
material parameters, independent of length. 



B. Electrostatically Dominated Limit 

The electrostatically dominated limit is defined by the 
condition 27rA2, 27rA3 3> Ai. Again, we investigate a tube 
under uniaxial stress only (/3 t = 0). In Fig. [3] we plot the 
equilibrium strain fields for several chiral angles. Note 
that in this case the tube is shorter than the strain decay 
length, and the strain fails to relax to its elastic limit. 
The tube develops significant torsional strain along its 
entire length unless it is in one of the high symmetry 
configurations. 

The bound charge density and the electrostatic poten- 
tial in the electrostatic limit are plotted in Fig. 2] for 
several chiral angles. The bound charge density is nearly 
linear throughout the interior of the tube and vanishes 
at the middle of the tube by symmetry. Near the bound- 
aries the charge density begins to develop a logarithmic 
divergence which is cut off because of the soft kernel. Not 
shown are delta functions in the charge density right at 
the tube boundary. Again, these are surface charge den- 
sities that arise due to our use of the soft kernel. The 
potential is also nearly linear along the entire length of 
the tube and vanishes at the center due to symmetry. 
Near the boundaries there are small, non-divergent cor- 
rections to the linear behavior. A linear electrostatic po- 
tential profile is what we would naively expect from the 
conventional 3-d case, but note the remarkably different 
conditions which give rise to the linear potential in this 
case. In the 3-d case the linear potential is produced by 
true surface charges interacting via a long-range electro- 
static interaction. In this 1-d case, the linear potential is 
produced by a nearly linear bound charge density inter- 
acting via a short-range electrostatic interaction. 

Inspection of the strains in Fig. [3] show they may be 
well fit by quadratic curves. An approximate analytic 
solution may be obtained in the form g s {x) ~ c Q + C2{x — 
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FIG. 3: Equilibrium solutions for a piezoelectric tube in the 
electrostatically dominated limit using the soft kernel. The 
upper panel shows the induced uniaxial strain for tubes of var- 
ious chirality, while the lower panel shows the corresponding 
induced torsional strains. The graphs were generated using 
(3 = 1, Ai = 1000, A 2 = 10000/tt, A 3 = 2A 2 ,and e 2 = v / 5000/ 7r " 
(charge/length). 
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This expansion is accurate to within a few percent pro- 
vided 7,r 3> Ai, which is generally true if A2, A3 1. 
The exception occurs if — > 7r/6, in which case t — > for 
any A2, A3. Even in this case the expansion differs from 
the numerical results only near the tube ends, and still 
gives acceptable agreement with the numerical results for 
the potential and total twist. For the potential across the 
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FIG. 4: Equilibrium solutions for a piezoelectric tube in the 
electrostatically dominated limit using the soft kernel. The 
upper panel shows the induced bound charge density for tubes 
of various chirality, while the lower panel shows the corre- 
sponding electrostatic potential along the tube. The graphs 
were generated using /3 = 1, Ai = 1000, A2 = 10000/-7T, 
A3 = 2A 2 ,and e 2 = V5000/7T (charge/length). 



tube we find 



AU = {(i s - c ) Ai 



12 Al 



(27) 



Because ((3 S — Co) 3> c-i the potential scales nearly linearly 
with tube length, an effect anticipated from our numeri- 
cal results. The total twist along the nanotube will also 
have a nearly linear scaling with length, provided the 
system remains in the electrostatic regime. 



IV. 



BORON NITRIDE NANOTUBES 



We now apply our theory to calculate the piezoelectric 
response of a BNNT under uniform longitudinal tensile 
stress. The 2-D elastic constants of a BN sheet are given 
by££ C11 .30 TPa nm and C 66 w .10 TPa nm, and 
the piezoelectric constant is given by£ e 2 = .12 e/Bohr. 
Experimentally accessible forces are on the nano-Newton 
scale; we arbitrarily assume a force of 1 nN. Typical tube 
dimensions are R — 1 nm and L — 1 /im. 

Using these values, we find that Ai = 1000, A2 = 
4.4 x 10- 3 , A 3 = 1.3 x 10~ 2 , /3 a = 5.3 x 10~ 4 and p t = 0. 
These numbers show that BNNTs are an extremely elas- 
tically dominated system, and our perturbation formulas 
given in Sec. (|III A[) apply. Equation (|2"Tj) agrees with 
the numerical results for the strain to better than 0.1%, 
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and Eqs. 1(22 ]) and (pS J) . with C v = 47re 2 sin(36»)/3 s /a and 
= 47rA 3 sin(3#) cos(39)(3 s /a, give the electrostatic po- 
tential and total twist with even better agreement. 

For all chiral angles the strain fields decay to within 
0.5% of their elastic limits at a depth of 2R, and are 
always within 3% of the elastic limit. The interior of 
the tube develops no significant twist and essentially re- 
sponds as an achiral nanotube. The total twist, in radi- 
ans, induced in the system as a function of chiral angle is 
A(j)(6) = 8.7 x 1(T 5 sin(30) cos(36>). The system develops 
significant surface charges and a bound charge density 
that decays rapidly on the interior of the tube, and the 
entire potential drop occurs within a few radii of the tube 
ends. Using the values from above and converting to SI 
units gives the potential across the tube as a function of 
the chiral angle as U{9) — 23.0sin(36') mV. 

Based on these calculations, BNNTs should develop 
a rather weak piezoelectric response with almost no in- 
duced twist regardless of chirality. Physically, this is ex- 
plained by the exceptionally large elastic constants of the 
BN sheet, which are inherited by the nanotube and make 
it a stiff system. BNNTs are so elastically dominated that 
the predicted strain decay lengths, between 1 and 2 tube 
radii, place this system near the limits of what can be 
described by our continuum theory. 



V. CONCLUSION 

In this paper we have developed and solved a contin- 
uum theory to describe the piezoelectric effect in 1-D 
tubular systems. These systems exhibit different physi- 
cal properties than their bulk counterparts. Unconven- 
tional physical effects include the development a bound 
charge density near the tube ends and a non-linear elec- 
trostatic potential profile along the tube. Furthermore, 
the piezoelectric response of these systems couples longi- 
tudinal strain and torsional strain. This coupling repre- 
sents a new degree of freedom which may be important 
in nanotube characterization and nanoscale device appli- 
cations. 

It would be interesting to measure the non-linear po- 
tential variation along a nanotube to verify the unconven- 
tional piezoelectric effects described here. BNNTs are the 
model physical system in which to realize these effects. 
However, BNNTs tend to be stiff structures for which 
the electrostatic contributions to the elastic properties 
are expected to be small. 

This work was supported by the Department of Energy 
under grant DE-FG-84ER45118. 
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